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The accelerated expansion of the universe has now been confirmed by several independent obser- 
vations including those of high redshift type la supernovae, and the cosmic microwave background 
combined with the large scale structure of the Universe. Another way of presenting this kinematic 
property of the Universe is to postulate the existence of a new and exotic entity, with negative 
pressure, the dark energy (DE). In spite of observationally well established, no single theoretical 
model provides an entirely compelling framework within which cosmic acceleration or DE can be 
understood. At present all existing observational data are in agreement with the simplest possibility 
of the cosmological constant to be a candidate for DE. This case is internally self-consistent and 
non-contradictory. The extreme smallness of the cosmological constant expressed in either Planck, 
or even atomic units means only that its origin is not related to strong, electromagnetic and weak 
interactions. Although in this case DE reduces to only a single fundamental constant we still have 
no derivation from any underlying quantum field theory for its small value. From the principles of 
quantum cosmologies, it is possible to obtain the reason for an inverse-square law for the cosmo- 
logical constant with no conflict with observations. Despite the fact that this general expression is 
well known, in this work we introduce families of analytical solutions for the scale factor different 
from the current literature. The knowledge of the scale factor behavior might shed some light on 
these questions since the entire evolution of a homogeneous isotropic Universe is contained in the 
scale factor. We use different parameters for these solutions and with these parameters we stablish 
a connection with the equation of state for different dark energy scenarios. 
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I. INTRODUCTION 



The research about the cosmological scale factor has attracted intense attention during the last years 
as well as the investigation about the connection between the cosmological constant and the dark energy 
issue. As the vacuum has a non-trivial role in the early Universe, a A-term in the Einstein field equations 
is generated. This A-term leads to the inflationary phase Based on the inflationary cosmology we can 
say that during an early exponential phase, the vacuum energy was a large cosmological constant. However, 
the current observed small value of the cosmological constant makes us to assume that A, representing the 
energy density of the vacuum, is a variable dynamic degree of freedom which being initially very large went 
down to its small present value in an expanding Universe f3|. The last one can be measured through the 
discrepancies between the infinitesimal value that the cosmological constant has for the present universe (it 
is very small in Planck units) and the values expected by the Standard Model 0]. 

There is an extensive literature that show examples of phenomenological A-decay laws. To mention some 
ideas, the studies comprise a A depending on temperature, time, Hubble parameter and scale factor 0]. A 
dynamically decaying cosmological constant with cosmic expansion has been considered by several authors 

ii,S0,ii[iE[iil,[i3- 

The importance of the subject resides in the fact that any nonzero value of A introduces a length scale and a 
time scale into the theory of general relativity The cosmological constant perturbs spacetime dynamics, 
although the general relativity works on scales much larger than the Planck scales. However, vacuum 
quantum fluctuations of the Standard Model rekindle the theory of the cosmological constant problem. 
Moreover, the cosmological constant is one of the candidates to be a DE model, namely, with the equation 
of state Lo = —I. A plaudible DE model is very important to explain the observations introduced by 



Supernovac (SNe) Type la Other candidates are scalar fields models such as quintessence, K-essence, 



tachion, phantom and etc.. For a review see [15| . 

Albeit there are many references about a scale factor inverse-square law for the cosmological constant there 
is a lack of information about the explicit form of the function of the scale factor that obeys this requirement. 



In reference there are some analytical solutions. However, those were obtained for a perfect-fluid matter 
energy-momentum tensor. Besides, it is our opinion that there is not an appropriate analysis in the literature 
that includes the scale factor analytical solutions in a DE approach. The underlying motivation to carry 
out these solutions is the fact that the time variation of A can lead to a creation of matter (i.e., A(t)CDM 
models) and/or radiation such as to help us in the investigation of the age of the Universe. 

In this work we will use the requirement of the energy-momentum tensor conservation in order to introduce 
some new analytical solutions for the scale factor different from the published and unpublished literature. 
While the current observations clearly pinpoint to a flat Universe, in the name of generalization, we will flnd 
solutions that have different curvatures, k = 0, —1, +1 as parameters, i.e., for flat, open and closed universes 
respectively. The obvious interest in a flat {k = 0) cosmological constant model has its motivation rooted in 
the fact that a A-term helps to connect inflation with observations. Moreover, with A, it is possible to obtain, 
for flat Universes, a theoretical age in the observed range, even for a high value of the Hubble parameter 



17|, [18|, but it is not the objective of this work. With these parameters we analyze some DE models in the 
background of these analytical calculations. In this paper we also assume, as will be said below, that the 
parameter of the DE equation of state, lu, is constant. This parame4ter measures the ratio between the dark 
energy pressure and its energyh density and is its precse value is one of the main tarets in general relativistic 
cosmology. 

In the next section we will depict the held equations and the differential form of the scale factor equation. In 
section HI we will present some motivations about the reasons to study this specific form for the cosmological 
constant. We present some possible analytical solutions for the scale equation in section IV. In section V we 
carry out a very brief review about DE and we analyze its connection with the obtained results. The final 
discussions are left for the final section. 
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II. DARK ENERGY 



The observations showed that the Universe was not only in expansion but that this expansion is not 
uniform, as beheved by the scientific community, but it is accelerated, One evidence of this acceleration is 
through the Hubble-Sandage diagram that shows the observed brightness of the SNe type la as a function 
of the redshift. It suggest that the bulk of the energy density in the Universe is repulsive and appears like 
an exotic form of energy with negative pressure. This exotic energy was dubbed as dak energy and it is 
considered as an additional component to the ordinary dark energy. There are many good reviews in the 



literature [15 1 



The analysis of dark energy phenomena begun with the observations of Supernovae (SNe) type la published 



independently by two different groups [14| in the late 90 's. As we said before the equation of state paramete 
to, is one of the most wanted number in relativistic cosmology today. One motivation for this quest is that 
if one could fix lo to be precisely —1, then, arguably, we could identify the DE with the vacuum state of all 
existing fields in the Universe, namely, the cosmological constant. On the other hand, if the precise value 
obtained is w ^ — 1, then we could not rule out the cosmological constant, but we can work on a model of 
the DE pressure responsible for the current cosmic acceleration as the potential energu density associated 
with a dynamical scalar field. In other words, we have the so-called quintessence field if —1 < uj < 1 
[? ] or a phantom component for w < —1 [35j. In both cases we have a violation of the strong energy 
condition {p + 3p > 0). In the case of the phantom component, moreover, we have the violation of the null 
energy condition (p + p > 0, for a homogeneous and isotropic Universe). Bounds on the value of uj can 
be obtained from observation such as the mentioned magnetic measurements of SNe la. There are other 
evidences that came with the anisotropics of the Cosmic Microwave Background (CMB) observations [3 



the Large-Scale Structure (LSS) data 31], the age of high-z objects [36|, measurements of the angular size of 



radio sources [37|, observation of the X-ray luminosity of galaxy clusters [38|, gravitational lensing [39] and 
there are others that the interested reader can find in [14] showing that a standard cold matter scenario is 
not enough. Recently, in [sit the authors showed, using the lookback time techique, how to provide bounds 
on the cosmological parameters with high accuracy. 

In spite of the number of DE models, the cosmological constant and the quintessence model are the most 
probables candidates. For the quintessence model, the equation of state is determined dynamically by the 
scalar field or tachyon. 

It is important to notice that the concept of an accelerated Universe and consequently of the existence of 
a DE component can be considered as an introduction of a new era in general relativity (GR) . It is possibloe 
that we have to review all the GR concepts since repulsive gravity is possible in GR. Through this section we 
will present briefly some possible alternative explanations for the acceleration of the Universe and connect 
these models with the results obtained in the last section, which is our actual objective. 

To accomplish this we have to analyze the main dark energy ingredient which is its equation of state 
given by equation p = lo p. We will see that the alternative models of dark energy have their origin through 
different values oi uj. It is worthwhile to stress that it is very plausible to consider this parameter constant 
because observational data can hardly distinguish between a varying and a constant equation of state [j^ . 
But there are many studies that consider this parameter as a function of the redshift z or the scale factor 
a{t). Here we will consider this parameter as constant. 

From the conservation of the energy momentum tensor through the Bianchi identities we have a continuity 
equation given by 

p = -3- {p+p) , (1) 
a 

and substituting the equation of state above and solving for p we have that 

p = p,a-^(^+-\ (2) 
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where po is the density of non-rclativistic matter today. For uj = we have that p = and this scenario is 
characterized by a dust dominated Universe, i.e., p — po a^^. For uj — ^ we have p — po a '^, i-e., a radiation 
dominated Universe. However, for both cases we have a deceleration expansion of the Universe scenario. 
Notice that it is the same result provided by the Newton gravitational theory. 

As we claimed in the first sections, the scale factor been proportional to the Universe radius, as said above, 
describes the growth of the Universe and it is defined so that at the present time to, a(to) = 1. It can be 
proved both classically and relativistically that the time evolution of the expansion obeys 

^=-^(p + 3,), (3) 
a 6 

where d{t) measures de acceleration of the Universe. 

The experiments described above confirm that at the present time, d(t) > 0, characterizing an accelerating 
Universe. This acceleration implies that, either: the Universe is dominated by some particles (dark matter) 
or by a field (dark energy) that has a negative pressure. In particular, substituting the equation of state in Q 
we have the condition u) < And again we have two situations: or there is in fact a non-zero cosmological 
constant; or the GR or the standard model, is incorrect. For w = — ^ we have an acceleration/deceleration 
boundary. The DE is so enigmatic that we can have two different scenarios: if DE is a cosmological constant 
the Universe will be acceleration indefinitely and will become empty; or if it is a dynamical scalar field, this 
last may decay and refill the Universe with matter and energy. 

Another scenario that appears in the literature is the one of a X-matter plus cold dark matter (the well 
known XCDM parametrization) where both fiuid components are separately conserved. This cosmological 
model is characterized by an equation of state — 1 < uj < 0. The model dominated by the cosmological 
constant p — — p, i.e., w = — 1 is a limiting case of the X-matter parametrization. In terms of energy density 
we have that, using p — pQ, and the energy density is constant in any epoch, hence, we can say that the 
cosmological constant corresponds to a fluid with a constant equation of state uj — ~1. 



III. THE SCALE FACTOR EQUATION 

Let us begin reviewing some basic concepts. The Einstein field equations are 

G^. + AffM- = SnGTf,, , (4) 
where A has often been treated as a constant of nature and where Gp,y is the Einstein tensor given by 

Gf_ii^ ~ R^u ^ '^RQtiv (5) 

and is the energy- momentum tensor. We are using relativistic units (i.e., c = 1) and we are assuming 
that the Universe is homogeneous and isotropic. In this case we will use the Robertson- Walker line element. 
Let us consider also a perfect- fiuid- like ordinary matter with pressure p and energy density p Q . 

Taking the covariant divergence of ^ and using Bianchi identities to guarantee the vanishing of the 
covariant divergence of the Einstein tensor, it follows (assuming the energy-momentum conservation law) 
that the covariant divergence of A g^i, must vanish also and therefore A is constant constituting a geometrical 
interpretation of A. 

As usual in recent works, we can move the cosmological term to the right side of Q in order to reconstruct 
the Einstein equation with an alternative form such as, 

G^, = 8^GV (6) 
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where 



A 
8^ 



(7) 



and we can interpret the cosmological term as part of the matter content of the universe, rather than a 
geometrical term. With this effective momentum-tensor T^i, satisfying the energy conservation as 

V f^, = (8) 

there are no reasons in order to refuse a varying cosmological terms At this point we have to observe 
that since we use c = 1, we will not consider here cosmological models with a time varying speed of light 

The incorporation of the cosmological term in the definition written in fT]) implies that the effective energy- 
momentum tensor has the perfect fluid form, with an effective pressure given by p = p — A/SttG and an 
energy density as p = p + A/8nG [i*]. 

Using the field equations (O and ([7]) and the law of energy-momentum conservation, equation ([5]), we can 
write 



SttG 



A 



pa + — a 

3 ^ 3 



and 



d_ 

da 



-3 p - 



A 

8^ 



(9) 



(10) 



For the equation of state we can write that 



P 



(7 - 1)P 



(11) 



with 7 = constant. Note from pTjl that for 7 = 1 we have matter-dominated epoch {p = 0) and for 7 = | 
we have radiation-dominated epoch. 

Notice that, in order to obtain a convenient form for athe nest equations we constructed the equation of 
state in pTjl . In a DE scenario, as said in the last section, we are clearly making uj — -f ~ 1. The reason for 
this will be clarified in the next lines. 

Substituting pT|) in pT?]) we have that 



Now, using ^ and (fT^ we obtain that 



1^^ dA 



SttG da 



37 
2 



A 



pa 



Using ([9]) to eliminate p in (fT3|) (or vice-versa) we can write the equation 



(12) 



(13) 



(14) 



which is a differential equation for the scale factor and has a cosmological term. We clearly see that this 
equation is independent of the fact that A is constant or not 

Our objective from now on is to present analytical solutions for the scale factor in equation But to 

accomplish this we have to choose a relevant expression for the cosmological term. 
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As said above, there is a whole hterature that approaches this issue (see [J] and the references therein). 
Here we choose a specific one (well known in the literature) but, we produce a variation in the curvature 
also. 

We can summarize the forms of the cosmological constant depending on the scale factor in a general 
expression as [g^] 

• \ 2 

a 



A{t) = aa-"' + I3[-] , (15) 



where a, /3 and m are constants. Ozer and Taha [6] and Chen and Wu [51 concluded, although in different 
contexts, that the time dependence of A should have /3 = and m — 2 20]. The interested reader can see 
the different variations concerning these parameters in 
Hence, let us choose the cosmological term as 

A = , (16) 

where B (we use here the same parameter used in 0]) is a pure number of order 1 which should be calculable 
in a model for the time variation of A. We can consider that the parameter ,B is a new cosmological parameter 
to be computed from observations, substituting A. If S = we have, obviously, an ansatz of A = 0. 



IV. SOME CONSIDERATIONS AND MOTIVATIONS 



With a time variation like the exposed in p5|) the values of A in the early Universe could be much bigger 
than the present value of A. However, even that order of small magnitude it was arguably large enough to 
drive various symmetry breakings that occurred in the early Universe. With a A = ansatz we cannot have 
such scenario. Otherwise, if ;B ^ 0, the equation (fT5|) can be identified as a "medium" time variation 

For a positive B our Universe is flat. In fact, considerations based on the second law of thermodynamics 
lead to the restriction B >0 6]. Also, if S > fc the Universe will expand forever; if S < fc the Universe will 
collapse in the future. For 7 > 1 in equation (|14p . even a closed Universe will expand forever and the k ^ 
case is no longer critical for the collapse of our Universe It can be shown that if our Universe is flat, 
observationally B should be positive [5|- 

If we write equation (|16p in a general form, i.e., as a function of the scale factor a~", like in the first term 
in equation (jl5p . the value n = 2 has dimensional reasons. With n = 2, both h and G disappear and we 
have, with c = 1 only, the equation (I16p [H]- The presence of H on the Einstein equations would be rather 
disturbing. 

Nevertheless, equation ()16p can be obtained theoretically based on general assumptions on quantum grav- 
ity. It can be considered also that this formulation has no conflicts with current data and can furnish some 
explanation about the inflationary scenario In other words we can say that this model is singular and 
preserves the standard picture of the early Universe. Notice that the condition dA/da < requires that 
B < independently of the curvature index k and hence A > for alH > [2^ . We will not make inflation 
considerations for the time being. 

In another point of view, the dependence on a used in ()16p was flrst introduced by Gasperini ^8:] in a 
thermal approach. An alternative motivation for this kind of dependence is that it also appear in some string 
dominated cosmologies [27|. We can also obtain an equation ([T6l) form based on semiclassical Lorentzian 



analysis of quantum tunneling 28[. Although this kind of dependence is only phenomenological and does 
not come from particle physics first principles J/7j] it can be considered as a laboratory in order to investigate 
the models it generalizes. As mathematically simple, we can obtain some analytical solutions. 
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V. THE ANALYTICAL SOLUTIONS 



After this brief historical review, let us now substitute in ([H)) to have that 



-.(1-2I)„.,.)J,,_2,,, ,17) 



which can be written as in Q as 



aa + Ai + A2 = , (18) 

where 

A. = -(l-|^)t^-|E, (19) 

and 7, fc and B constitute a group of free parameters which have the respective physical considerations 
mentioned above. In our analysis, the parameter fc, as said above, will have the values 0, —1 and +1. Just 
to remember, if we observe (IT51) we can see clearly that if S = 0, then the expression covers the A = case. 
If K ^ [1] we have two options, B > k where the Universe will expand forever ot B < k and then the 
Universe will collapse in the future. If S > 1, even a closed Universe (with k — 1) will expand forever and 
the fc = case is no longer critical for the collapse of our Universe [1]. 
For the case Ai ~ —1/2 and substituting in we have that 

(Pa \ f da 
~ 2 V 



^^--A-] +A2 = 0. (20) 



After some algebraic work we can reduce the solution to a simple form as 

a{t) = - 2 ^ + ^ (t + C2f , (21) 
ci 4 

where Ci and C2 are constants to be determined. Observe that we have written the equation above in order 
to maintain the parameter A2 in evidence for the following analysis below. 

For this case of Ai = —1/2, from equations ([TO]) we have that 7 = 1/3 and consequently that A2 has a 
dependence on k and B that can be written as 



K2 = --k --B , (22) 
2 6 ^ ' 



1,-1, 






Ar° = - 


> 






1 

2 " 




Ar^ = - 


1 

" 2 " 





(23) 



and the analysis about the expansion or collapse of the Universe can be accomplished as above. 
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Since the value B = can be obviously discarded, we can stress that the general behavior of the solution 
(PT|) does not change. As a final observation we can say for Ai = —1/2 the equation of state can be written 
as p = — 2/3, where we have an inflation scenario. 

As usual we measured in units of Hubble times (r = Hot) and we have imposed the following boundary 
conditions at present epoch: 



a (to) = 1 and H (tq) = Hq. 



After have calculated the expression for the Hubble parameter we used these boundary conditions to 
determine the free parameters values in the scale factor since we have two equations (scale factor and 
Hubble parameter) and two incognites cl and c2. 

To plot the figures of this work we have set Hq the 74 that is inside of the current accept value of Hq (73 
±10 Km s Mpc 1) 0. 

In figure 1, after the change of variable r = Hot, we present the plot of the scale factor versus r with 
Ai = i and A2 varying from —20 to 20. 




FIG. 1: Scale factor, equation (18), versus r with Ai = | and A2 varying from —20 to 20 (Units Hq ^ 



For Ai = — 1 we have two general solutions for the equation (fT7|) . 



1 

ai{t) - — 
4ci 

1 

4^ 



a2{t) 



gCl(t + C2) —4^9 6"^^'* + '^^' 



(24) 



where the difference is the signs of the exponentials. 

For Ai = — 1 we have a little different behavior. From (fT9|) we have that 7 = and consequently that 



A, 



-k, so we can write that A2 — 0, A2 



k=l 



4 and A 



k=-l 



= 1. For A2 = the solutions (l24l) have 



different behaviors altogether, namely. 



1 

ai(t) = — e 

Cl 

1 

02 (i) = — e 

Cl 



Cl(t+C2) 
Cl(t + C2) 



(25) 



and the equation of state now is p = 



—p, again an inflation scenario. 



For A2 = ±1 we have that, 



,fe=±i 



,fe=±i 



,Cl{t + C2) 



-Cl(t+C2) 



it) 



4ci 



g-Cl(t+C2) -p 4gCl(t+C2) 



and we can see again completely different formulations for the scale factor in an inflation scenario. 



In figure 2, we present the plot of the scale factor versus r, first equation (23), with Ai = — 1 and A2 
(fc = 1 and 7 = 0) . 




FIG. 2: Scale factor versus t, first equation (23), with Ai = 



— 1 and A2 = 



-1 (fc = 1 and 7 = 0) (Units Hq 
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In figure 3, we present the plot of the scale factor versus t, second equation 23, with Ai = — 1 and K2 — —\ 
(fc = 1 and 7 = 0) . 




FIG. 3: Scale factor versus r, second equation (23), with Ai = — 1 and A2 = — 1 (fe = 1 and 7 = 0) (Units ffg ^ ) 

In the figure 4, we present the plot of the scale factor versus r, first equation (23), with Ai = — 1 and 
A2 = 1 (fc = -1 and 7 = 0) with singularity in r=-0.1284233275. 




FIG. 4: Scale factor versus r, first equation (23), with Ai = — 1 and A2 = 1 (fc = —1 and 7 = 0) with singularity in 
r=-0. 1284233275. (Units H,^ ) 



For Ai = 1 we have the following solution. 



a{t) = ±J|L _ A2(t + c,Y , 



(27) 



and from equations we have that 7 = 4/3 and 



^2 = k- -B, 



(28) 



and the general behavior of the solutions (P7)) do not change since we have the condition B ^ Q. However, 
we have to notice that the square root demands attention in order to not obtain imaginary scale factors for 
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fc = 0, ±1. But this is a trivial task to substitute in the values of A2 obtained in (pS)) with different k. 
We will not have any new physical feature in this specific mathematical analysis. The equation of state is 
p = ^ p, a radiation-dominated epoch as said above. 

Next we show a solution that does not have constraints in Ai and A2 besides those given by the equations 
(16) 

ait) = -^^""^(^^+^^^»'. (29) 

If we set Ai = —1 and A2 = 1 the solution has a singularity in r = —0.1279126267. Figure 5 show this 
result. 




-0.15 -0.1 -0.05 

Time tau 



FIG. 5; Scale factor versus r, equation (26), with Ai = — 1 and A2 = 1 (fc = — 1 and 7 = 0) with singularity in 
r=-0. 1284233275. (Units ) 



VI. FINAL CONSIDERATIONS 



The equation was studied in 

0, S, @, mi, 0, Hi, and it can also be obtained theoretically fr 



rom 



some assumptions from quantum cosmology. It is not in conflict with experimental observations and with 
the inflationary scenario. From quantum cosmology arguments it is more convenient to use the scale factor 
instead of the age t of the Universe Q . 

In this paper we have computed some analytical solutions for the scale factor since the effective cosmological 
constant A is time varying according to an inverse square law in the scale factor. This inverse square law can 
be justified based on the principles of quantum cosmology. Since such time variation of A leads to creation 
of matter and/or radiation, we believe that the formulation of a scale factor analysis is of current interest 
both theoretical and experimental. Another motivation can be found in the strict relation between the 
cosmological constant problem and the problem of the age of the Universe treated in the light of quantum 
cosmology. 

There is a great literature about different functions for the cosmological constant. Specifically concerning 
a function like in the equation (jl6p the existing published papers do not bring any analytical expressions for 
the scale factor itself as a building term in A. 
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As we are interested in a general analytical behavior of A, numerical solutions of psp are out of the scope 
of this work. Therefore we found three analytical solutions. Other alternatives varying the parameter Ai in 
(|18p brought singularities which demands numerical efforts. 

Our main motivation is to show simply that there is another possible and different approach concerning 
the cosmological constant in order to obtain some progress in some questions about the Universe. Although 
the model proposed by Chen and Wu in equation (|16p are compatible with the astrophysical data, it alters 
the predictions of the standard model for the matter-dominated epoch. 

As a future perspective it can be investigated, with the solutions obtained here, some aspects of inflation 
like the current value of the deceleration parameter and the parameter Oqj namely, to obtain its value and 
test if they are compatible with the standard model demands. Note that the knowledge of the form of the 
scale factor, as well known, can describe if an Universe is an expanding one or not through the line element. 

An investigation about the age of the Universe in different scenarios as well as the equations of state, with 
the results obtained in this paper, is a work in progress. 
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Appendix 



Full List of Scale Factor Solutions. 



Quadratic solution: 



a{t) = -2^ + 2l^t + c,)\ 
ci 4 



Exponential type solutions: 



ai{t) 



4ci 



If Ai = 1, with square root type solution: 



gCi(t+c2) _ 4 A2 e"'^^^*"'"'^^' 



a{t) = ± - A2(i + C2)2 , 



Solutions with arbitrary Ai and A2: 



a (t) = C3 ± 



V-A1A2 (ci +C2t) 
A1C2 



If Ai = i (n e iV), where n is a positive even number, the solution is 

a (i) = y^A^A7(t+^ 



If Ai = — 1, we have the following solutions 



a{t)=± 



7^ cosh {C2 + C3t) 

C3 



u\ - 1 , ^/=A;(cosh(C2 + C3^)r 

"^*^"^2 C3 ^ C3 



, . _ V=A2COSh(C2 + C3t) , 4 V=A2(cosh(C2 + C3^)) 
"^^"^ C73 ^3 C3 



a(f) = ±- 



/A^sin ((72 + ^3^) 
C3 



1 VA^ , VA^(sin(C2 + C3^))' 



\/A^sin(_Cg + _(7gt) , 4 ^/A^(sin {.C2 + .C3t)y 

~ ^ IC^ 3 :c3 



a{t) = ±- 



ZC2 ' 



„ ^ - ^ ^ ^ + (-^^ + *))^ 



y^A^cos (_CJ + i) 4 \/=A^(cos (.Ci + .C2 t)f 
" ~ ^ :C2 3 leg ' 



a{t) = ± 



V=A2Cosh(_Ci + .C^t) 
IC2 ' 



i V^-^ . (cosh (-C^ + -C^ i))^ 



2 .C^ 



.C2 



cosh {_C1 + .C2t) , 4 ^/^{cosh{.Cl + .C2t)f 
" " ^ leg 3 ' 



a{t) = ±- 



ZC2 ' 



, . _ 1 VA^ y/A^ {sin {_C1 + .C2t)y 



x/A;sin(_Ci +_C2t) , 4 VA;(sin(_CT 
''^*>-^ _C2 3 .C^ 



a(t) = ±- 
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y/M sinh {_C1 + .C2t) 4 y/M{smh {_C1 + .C2t) f 

:c2 ^3 :c2 ■ (2^) 



If Ai = —2, wc found out Jacobi elliptic type solutions 



1 ^/TX^JacohiDN (_C2 + .C3 t,y/2) 
a{t) = ±- L, (24) 



^f^JacohiNC (_C2 + .C3 1, i^/2) 
a{t)^±^ ^ (25) 



1 ^f=2K^JacobiND LC2 + .C3t,V^ 
ait)^±- -^^ (26) 



_ 1 ^/'^2MJacobtNS {.C2 + .C3 1, i) 
a(i)-±2 :C3 ' ^ ' 



«W-±2 • (28) 
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